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In this paper we study an s-wave topological superconductor (SC) with a square vortex-lattice. 
We proposed a topological Majorana lattice model to describe this topological state which was 
supported by the numerical calculations. We found that the Majorana lattice model is really a 
"topological SC" on the parent topological SC. Such hierarchy structure becomes a new holographic 
feature of the topological state. 



I. INTRODUCTION 

Recently, the search for exotic states supporting Majo- 
rana fermions (modes) has attracted increasing interests 
due to their potential applications in fault-tolerant quan- 
tum computations [lMl 5| . A creative proposal is the prox- 
imity effect between s-wave superconductor (SC) and 
topological insulator Q. This system exhibits non-trivial 
topological properties, including the nontrivial Chern 
number in the momentum space, the chiral Majorana 
edge states, in particular, the Majorana mode around 
the 7r-flux vortex. Another possible example of such 
quantum exotic states is the chiral p + ip topological 
superconductors Q. The quantized magnetic vortex (tt- 
flux) in the chiral p + ip topological SC hosts the Majo- 
rana zero modes and obeys non-Abelian statistics. When 
there are two 7r-fluxes nearby, the intervortex quantum 
tunneling effect occurs and the Majorana modes on two 
7r-fluxcs couple. The tunneling amplitude is determined 
by the overlap of the wave function of the localized Majo- 
rana bound states [3 [l^. Thus, such tunneling must be 
taken into account when the average distance between 
localized 7r-fluxes becomes the order of the Majorana 
bound state decay length. For the topological quantum 
computation based on the non-Abelian anyons, the tun- 
neling effect would split the zero-energy bound states and 
lift the ground state degeneracy. Beside, the sign of en- 
ergy splitting is important for understanding the many- 
body collective states [l2i [Trj Il7j|. 

However, till now people have not identify the exper- 
imental realizations of p-wave SC in condensed matter 
physics [§l [20| . Instead, people pointed out that the chi- 
ral p-wave SC may be realized in ultra-cold fermionic 
atom systems [2l| - |2a |. Theoretically, p-wave SC may be 
realized via a p-wave Feshbach resonance in experiment. 
Due to huge particle- loss [14|, this idea also has not been 
realized. For this reason, people proposed a more advan- 
tageous scenario based on s-wave SC of cold fermionic 
atoms with laser-field-generated effective spin-orbit (SO) 
interactions and a large Zeeman field Q. In this scenario, 
there is a non-Abelian topological phase that is different 
from p-wave SC, in which the SO interaction plays the 
role of the p-wave SC order parameter. 

In this paper, we start from this s-wave topological SC 
model and show that there exists the Majorana mode 
hosted by the 7r-fhix (quantized magnetic vortex). Then 



we focus on the super-lattice of 7r-fluxcs. Because each ir- 
flux traps a Majorana mode, we can have a lattice model 
of the Majorana modes that is called the Majorana lattice 
model. We took into account of zero mode tunneling that 
couples the vortex sites. We found that this model shows 
nontrivial topological properties, including a nonvanish- 
ing Chern number, chiral Majorana edge state [2(1] . In 
this sense, the Majorana lattice model is really a "topo- 
logical SC" on the parent topological SC. Such hierarchy 
relationship between the Majorana lattice model and the 
s-wave topological SC model with vortex-lattice is a new 
holographic feature of the topological states. 

The paper is organized as follow: We first introduce the 
topological s-wave SC model on a square lattice in Sec. 
II, then study the topological properties of this model. 
In Sec. Ill, the Majorana mode around a 7r-flux is ob- 
tained and the intervortex quantum tunneling effect is 
also studied. We also study the topological SC with a 
square vortex-lattice numerically. In Sec. IV we write 
down a Majorana lattice model to describe the coupling 
effect between Majorana modes trapped in vortices, and 
the topological properties of this Majorana lattice model 
are also analyzed. Finally we draw the conclusion in Sec. 
V. 



II. THE S-WAVE PAIRING TOPOLOGICAL 
SUPERCONDUCTOR WITH RASHBA 
SPIN-ORBITAL COUPLING 



As a starting point, the s-wave pairing SC with Rashba 
SO coupling is defined on a square lattice which is 
described by 
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where the kinetic term Hk, the Rashba spin-orbital (SO) 
coupling term H SOl and the superconducting pairing term 
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where c J(T (cj CT ) annihilates (creates) a fermion at site j = 
(jx,jy) with spin a = (t,i), /J- = x or y, which is a basic 
vector for the square lattice. A serves as the SO coupling 
constant and A as the SC pairing order parameter. In 
addition, in order to observe nontrivial phases supporting 
Majorana modes excitation, the chemical potential u and 
the Zeeman term h also be included in this model. 

The Hamiltonian T-L can be transformed into momen- 
tum space from cj = l/\/L'^2 k Cke~ lkR j . Writing Cj in 

the particle- hole basis ip\ = (ckt, Ck: ■ 
tain its Bogoliubov-de Gennes(BDG) form 
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where the Bloch Hamiltonian H(k) is a 4 x 4 matrix 
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e(k) = —2t s (cosk x + cosk y ) — u 
g(k) = — 2A(sinfcj, + isuik x ). 
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Then the energy spectrum e± is obtained by diagonaliz- 
ing H(k) as 



e± =±[m 2 (fc)+| 5 (fc)p 



-h 2 ±2y/\g(k)\^(k) 
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where m(k) = y / e 2 (k) + A 2 . 

In FigfTl we plot the global phase diagram. The blue 
lines are the boundaries between the topological SC and 
the SC with trivial topological properties. In the topo- 
logical SC, the ground state has a nontrivial topological 
invariant (Chern- number) Q = lQ. Such topological 
invariant is robust, for the Rashba SO coupling can be 
mapped into a p + ip SC gap through a suitable unitary 
transformation (see details in Ref.@). We also study the 
density of state (DOS) in the topological superconductor 
numerically and show the results in FigJ5J From this re- 
sult one can see that there exists a finite energy gap of 
the topological SC state. In other regions of the phase 
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FIG. 1: (Color online) The phase diagram (A = A = 0.5i s ): 
The blue lines are the boundaries between the Q = 1 topolog- 
ical SC (TSC) and the normal SC (NSC) with trivial topo- 
logical properties. The color denotes the size of the energy 
gap. 




FIG. 2: (Color online) The density of states (DOS) of the s- 
wave SC model for the case of u = — 4t 3 , A = 0.5t a , A = 0.5t s , 
and h — 0.8t 3 . There exists an energy gap AE = 0.7t„. 



diagram, the ground states are the SC state with trivial 
topological properties. 

The operators of the Majorana edge states of this topo- 
logical SC as shown in Fig. 3 satisfy 7q = 70 , which 
indicates the edge states are really Majorana fermions. 
Indeed, we can see that the particle-hole operator C = 
a x ® 1 acts on %{k) as 



CU(k)C- x = -U*{-k) 



(7) 



implies that the Bogoliubov quasi-particle operator fol- 
lows 7l fe = 7fc . From Fig. 3 (a), one can find that the 
Majorana edge modes cross zero energy at k y = for 
the topological SC with u = — 4t s , A = 0.5t a , h = 0.8t s . 
The odd number of crossings leads to the topological pro- 
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FIG. 3: (Color online) The spectral flow of Hamiltonian H on 
cylinder with u — — At s , X = 0.5f s , A = 0.5f s , and h = 0.8t s 
for (a) but h = 0.4f s for (b). For topological SC, there exist 
chiral gapless edge states in (a); while for normal SC there is 
no edge state in (b). 



FIG. 5: (Color online) The energy splitting SE as the func- 
tion of the spacial distance of two 7r-flux- vortices, decays ex- 
ponentially and shows the oscillating behavior with D. The 
parameters are chosen as it = — 4t s , X — 0.5t a , h = 0.8t a , 
A = 0.5t 3 . a is the lattice constant. 



III. MAJORANA ZERO MODES AROUND 
VORTICES OF THE TOPOLOGICAL SC 



A. Majorana zero modes around a pair of vortices 
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FIG. 4: (Color online) The particle density distribution of 
two localized zero mode states around the two 7r-flux-vortices 
for the case of u = — 4t 3 , X = 0.5t s , h = 0.8t 3 , A = 0.5t a . 



tected Majorana zero modes. 
u = -At s , X = 0.5t SJ h = QAt s 
(Fig.3 (b)). 



For the trivial SC with 
there is no such crossing 



We now start with the discussion on the Majorana zero 
modes which associate with a pair of 7r-flux-vortices. The 
particle density distribution of the electrons is shown in 
FiglU The main result shows that there appear two 
approximately zero modes in the presence of two well- 
separated 7r-flux-vorticcs. When two 7r-fiuxes are well 
separated, the quantum tunneling effect can be ignored 
and we have two quantum states with zero energy. On 
the other hand, for the small spatial distance D, the vor- 
tices interaction becomes stronger and the energy split- 
ting can not be neglected. As shown in FigJSJ the energy 
splitting 8E as the function of D, oscillates and decreases 
exponentially. Likewise, a pair of 7r-flux- vortices in the 
topological p-w&ve SC or quasi-hole in the Moore-Read 
state has the qualitative similar behavior fiH li~9l. Wf\ . 



B. Topological properties of the topological SC 
with a square vortex-lattice 

Next, we study the topological SC with a square 
vortex-lattice numerically. The illustration of the vortex- 
lattice (7r-flux-lattice) was shown in FigJBl 

The tunneling effect between vortices would lead to 
the coupling between two Majorana fermions around the 
vortices. Thus there exists a mid-gap energy band for 
the Majorana fermions. We study the density of state 
(DOS) in the topological SC with a square vortex-lattice 
numerically and the results are shown in FigjTJ From 
FigJTl one can see that besides the energy bands of the 
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FIG. 6: (Color online) The illustration of vortex-lattice: the 
distance between two nearest vortices is 6a. The parameters 
are chosen as u = — 4i s , A = 0.5t s , h — 0.8t s , A = 0.5t a . 
Each vortex traps a Majorana mode. Thus, in Sec. IV we 
have a tight-binding Majorana lattice model to describe the 
low energy physics of the multi-Majorana-mode. 




FIG. 7: (Color online) The DOS of topological SC with square 
vortex-lattice: the mid-gap states are induced by the vortex- 
lattice. The inset shows the details of the mid-gap energy 
band. 



paired electrons there exists a mid-gap energy band in 
parent topological SC. In particular, the mid-gap energy 
band has finite energy gap. It means that this mid-gap 
system as shown in Fig [7] may be a topological state asso- 
ciated with a non-trivial topological number intuitively. 
To check the topological properties of the mid-gap sys- 
tem, we calculate its edge states. We consider a system 
on a cylinder with 12 super-unitcells along x-direction 
while periodic boundary along y-direction. Thus, k y is 
still a good quantum number and permits the Fourier 



transformation Ck y (j x ) 



E 7 „ c(j x ,j y )e lk yiy. The 



/L y 

spectral flow of this system on a cylinder is plotted in 
FigJHl From FigJSl one can see that two gaplcss chiral 



FIG. 8: (Color online) The spectrum flow of topological SC 
with square vortices lattice on a cylinder. There exist gapless 
edge states for the case of u = — 4t s , h = 0.8t s , A = 0.5i s , 
A = 0.5i s . 
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FIG. 9: (Color online) The spectrum flow of normal SC with 
square vortex-lattice on a cylinder. The parameters are u = 
— 4t s , h = 0.2t s , A = 0.5t 3 , A = 0.5t s . There is no edge state. 



edge states are localized at the boundaries. On the other 
hand, we also plot the spectral flow in FigO When the 
parent SC is a non-topological SC, the edge states disap- 
pear. As a result, we conclude that the mid-gap system 
shown in Fig{7] is indeed a topological state. From Fig. 9, 
we also observe the disappear of the mid-gap band which 
is induced by the vortex-lattice. So the Majorana mode 
around the 7r-flux is protected by topological invariant of 
the parent s-wave SC. 



IV. THE TOPOLOGICAL MAJORANA 
LATTICE MODEL 

In the last section we found a topological mid-gap sys- 
tem of the s-wave topological SC with a square vortex- 
lattice numerically. In this section we will learn the na- 



ture of this topological mid-gap system analytically. 

We propose an effective description of the s-wave topo- 
logical SC with a square vortex-lattice, of which each vor- 
tex traps a Majorana zero mode and two Majorana zero 
modes couple with each other by a short range interac- 
tion. The interaction strength is just the energy splitting 
5E from the intervortcx quantum tunneling. We call this 
effective description as the Majorana lattice model, of 
which the tight-binding Hamiltonian can be written as 

Um.f = i s 3ktjklklj (8) 

where tjk is the hopping amplitude from j to k, and satis- 
fies tjk = t* k . 7j is a Majorana operator (7J = jj) obey- 
ing anti-commutatc relation {7j,7fc} = 28 jk- Sjj = — Sji 
is a gauge factor. The pair (j, k) denotes the summation 
that runs over all the nearest neighbor (NN) pairs (with 
hopping amplitude t) and all the next-nearest neighbor 
(NNN) pairs (with hoppin g am plitude if). From the poly- 
gon rule proposed in Ref. [28| . each triangular plaquette 
possesses it/2 quantum flux effectively. This Hamiltonian 
allows a Z2 gauge choice Sjk = ±1. For our Majorana lat- 
tice model, one possible gauge is given by FigJTDJ Thus, 
the total number of Majorana modes N must be even 
and then we can divide the Majorana lattice into two 
sublattices. 

Then the Majorana modes can be combined pairwise 
to create N/2 complex fermionic states by pairing the 
Majorana operators as 

7 2 y 1 = aj + at, j b 2j = {a 3 - a])/i (9) 

where a,j (at) annihilates (creates) a fermion at link j. 
In terms of operator aj, the Majorana Hamiltonian takes 
the form of a " topological SC" state as follows: 

n m j = ^ [*(4 a x+j - Oxfll+j) + t{a x a x+i - 4«4 +j ) 

X 

+ 2t(a x a x+i - a x a x+i ) - 2it(a x a x+i + o+a^ +i ) 

+ t'(al a x+i _j - Oxai+i-j - axflx+i-j + aiai + i_j) 

+ t'(al a x+i+j - a x a x+i+j + flxflx+i+j - ax a I+i+j) 

- t(a x a x - a x al)] (10) 

where i and j are two orthogonal unit vectors (FiglTOj) . 

Performing a Fourier transformation ak = 
—h= y\ a x e _4fc - Rx , the Hamiltonian becomes 

V L x ^— 'X 

Hm.f = J2lptd a (k)T a lP k (11) 

k,a 

in the basis ip^ = (al,a_fe)) where a = x,y,z, r a is the 
Pauli matrix, and 

d z (k) = — 2t sin 2 k y + At' cos k x cos 2 k y 

d y (k) = ~sm2k y (t + 2t'cosk x ) (12) 

<F{k) = ~2tsink x 
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FIG. 10: Two sublattices of the Majorana lattice: 7" and 
7 6 . The arrows denote the e" 1 ^ 2 of the hopping parameters. 
Accordingly, there exists n/2 in each triangular plaquette. 
The lattice constant of the Majorana lattice model is b — 6a. 
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FIG. 11: (Color online) The DOS of Majorana lattice model 
for the case of t = 0.048f s and t' = 0.006i s . 



See the detailed calculations in Appendix. Then we get 
the energy spectrum E{k) of the Majorana lattice model 
as 



E{k) = ± i^wi 2 - ( 13 ) 

From this result we can derive that there always exists 
an energy gap of the Majorana lattice model as long as 
t' 7^ 0. Then we calculate the DOS of the Majorana 
lattice model and show the result in FiglTT] One may 
compare the DOS of the Majorana lattice model in FigfTTI 
and the DOS of the mid-gap system in FigJ7]and find the 
similarity between them. 

To characterize the topological properties of the Ma- 
jorana lattice model, we define the nontrivial topological 
invariant - the Chern number 

which measures that the unit vector d(fc)/|d(fc)| maps 
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FIG. 12: (Color online) The edge state of Majorana lattice 
model for the case of t = 0.048i s , t' = 0.006i s . 

the Brillouin zone boundary onto sphere S 2 via the Chern 
number Q. In the presence of NNN hopping t', we have 
Q = ±1. Furthermore, we calculate the edge states of the 
Majorana lattice model and present the result in FigfT2] 
One may compare the spectrum of the edge states of the 
Majorana lattice model in FigJT2]with that of the mid- 
gap system in FigJ5]and also find the similarity between 
them. Now we can conclude that the (topological) Ma- 
jorana lattice model captures the key low energy physics 
of the topological SC with a square vortex-lattice. 

In addition, we finish this section by a brief discussion 
of the s-wave topological SC with triangle vortex-lattice. 
If in the region of \t'/t\ << 1, we observe that the model 
is characterized by the winding number Q = ±1, and 
this may stem from the redundant NN hopping (which 
behave as the NNN hopping of the square Majorana lat- 
tice model) of the triangular lattice. Moreover, the exotic 
phase Q = ±3 can be reached by varying of the ratio t' /t. 
In the region of \t' /t\ >> 1, the triangular Majorana lat- 
tice model has a 2 = ±3 phase. It's worth to point out 
that, for the triangular Majorana lattice model in p + ip 
SC, or interacting vortex-lattice in Kitaev's hone yco mb 
model, people may find similar phase diagram(26l l29j] . 



V. CONCLUSION 

In the end, we draw a conclusion. In this paper we 
studied the properties of an s-wave topological SC with 
a square vortex-lattice. Because each vortex traps a Ma- 
jorana zero mode, when we took into account of zero 
mode tunneling that couples the vortex sites, the Ma- 
jorana zero modes of the vortex-lattice form a Majorana 
lattice model. We found that this Majorana lattice model 
shows nontrivial topological properties, including a non- 
vanishing Chern number, chiral Majorana edge state. In 
this sense, the Majorana lattice model is really a "topo- 
logical SC" on the parent topological SC. And the "topo- 



logical SC" induced by the vortex lattice is topologically 
protected by the topological invariant of the parent s- 
wave SC. Such correspondence between the topological 
properties of the Majorana lattice model and the topolog- 
ical properties of the s-wave topological SC model with 
vortex-lattice is another holographic feature of the topo- 
logical state. In addition, we also had used the numeri- 
cal approach to study the s-wave topological SC with a 
square vortex-lattice and got similar results. 

Furthermore, the Majorana lattice model for chiral 
p + ip topological superconductor with vortex-lattice and 
that for the coupled system due to the proximity effect 
between s-wave SC and three dimensional topological 
insulator with vortex-lattice have the same topological 
properties to our case. Thus, this approach paves a new 
way to observe the Majorana modes for the topological 
SC. 
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Appendix A: Fourier transformation of the 
Hamiltonian for the Majorana lattice model 

The Hamiltonian of the Majorana lattice model 
can be obtained by Fourier transformation a x = 
l/V^E k ake~ 4ki S 

n m . f = i/L^Yl 

k,k' x 

{i(4a k< e ik ^- ik '^+J - a k 4,e- ikfl * +ik '^+J) 
+ t(a k a k ,e- ik ^- ik ' fl -+J - a k 4,e ik ^ +ik '^+J) 
+ 2i> k a k 'e ik ^- 4k ' fl -+ i - a^e"^^'^) 

- 2it(a k ave- lWR *- lk ' R ^ + a k a k ,e jki? * +ik ' 

+ t'(aia k ,e tkR *- lW ' R ^-> - a k al,e- ikR - +ik ' R -+^) 

- t'(a k a k ,e- ikR -- lk ' R -+--i - a+a^e^+^+'-J) 
+ t'ta^e^-^+'+J - a k a k ,e- ikii * +ik ' 

- t(ala kl e ikR '- ik ' R " - a k ale~ lkR - +tk ' R -)} 

(Al) 
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where i? x +<5 = Rx + S, (S = i,j, i + j). Using the identity 
<Sk,k' = l/ixEx ei(k " k ' )Jix . wc havc 

n m .f = J^iM^e-^ - a k ale ik *) 

k 

+ t(a k a- k e? ky - a\<J_ k e- ik *) 
+ 2t'(ala k e- ik * -a k ate ik *) 

- 2it{a k a_ k e lk * + ala i _ k e~' lk *) 

+ t'(ala k e- i(k *- k y^ - a k ale i{ - k "- k ^) 

- t' (a k a- k e i{k,c ~ ky) - a\a}_ k e-^ K - k ^) 
+ t'{ala k e~^ +k ^ - a^le 1 ^^) 

+ t'{a k a. k e i ^ +k ^ - a{al k e^ k * +k ^) 
-t{a\a k -a k a\)} (A2) 



which is 

Wm.f = (t cos ky + 2t' cos k x + 2t' cos k x cos ky)a\a k 

— (i cos fc y + 2t' cos fc^ + 2t' cos fc^ cos k y )a k a\ 

— (2t sin k x + it sin k y + 2it' sin fc y cos k x )ci-k<i k 

— (2t sin fc^ — it sin fcj, — 2it' sin k y cos 

(A3) 

or 

?W = E E ^ a (k)^a^k- (A4) 

k a—x,y,z 

where ^ = (aL a -k)- 
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